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From the microscopic quantum Langevin equations 
(QLEs) we derive the effective semiconductor QLEs and the 
associated noise correlations which are valid at a low-injection 
level and in real devices. Applying the semiconductor QLEs 
to semiconductor light-emitting devices (LEDs), we obtain a 
new formula for the Fano factor of photons which gives the 
photon-number statistics as a function of the pump statis- 
tics and several parameters of LEDs. Key ingredients are 
non-radiative processes, carrier-number dependence of the 
radiative and non-radiative lifetimes, and multimodeness of 
LEDs. The formula is applicable to the actual cases where 
the quantum efficiency n differs from the differential quantum 
efficiency rjd, whereas previous theories implicitly assumed 
Tj = rjd. It is also applicable to the cases when photons in 
each mode of the cavity are emitted and/or detected inho- 
mogeneously. When r/d < rj at a running point, in particu- 
lar, our formula predicts that even a Poissonian pump can 
produce sub-Poissonian light. This mechanism for genera- 
tion of sub-Poissonian light is completely different from those 
of previous theories, which assumed sub-Poissonian statistics 
for the current injected into the active layers of LEDs. Our 
results agree with recent experiments. We also discuss fre- 
quency dependence of the photon statistics. 

PACS number: 42.50Lc, 42.50Lv 



I. INTRODUCTION 

There have been many researches on quantum noise in 
light-emitting devices (LEDs) since the celebrated work 
of Shawllow and Townes [Q . We here consider the quan- 
tum noise in semiconductor LEDs. This subject was first 
studied by Haug and Haken [^) , who derived from a mi- 
croscopic model useful formulae for the first and second 
optical coherence of LEDs. 

Recently, much attention has been paid to sub- 
Poissonian light (SPL), which has the intensity fluctu- 
ations lower than the standard quantum limit [g| , of var- 
ious LEDs j|-|l4| • The mechanism for generation of SPL 
in LEDs is quite different from that of squeezed light 
in nonlinear crystals The latter mechanism is well 
understood as a Bogoliuvov transformation of a coher- 
ent state into a squeezed state of light jll| . The former 
mechanism, on the other hand, is often described by the 
Langevin theory of laser PJl7|-|2l|. Previous studies on 
SPL in LEDs assumed that only a single mode or a few 
modes of photons are excited j4||j,[l4| . When the injection 
level is lowered, however, many modes of photons become 



relevant, and we must consider all of them. Simple theo- 
ries for such a case were reported in Refs. HfjJ. However, 
they are too simplified so that they cannot explain recent 
experiments by Hirano, Kuga (HK) and coworkers [p[~pd], 
who demonstrated that the experimental results at a low- 
injection level (L1L) disagree with the predictions of the 
simplified theories. The disagreement appears when the 
quantum efficiency rj differs from the differential quan- 
tum efficiency 7/d. HK |J suggested that non-radiative 
processes might be responsible for the disagreement. 

In this paper, to resolve the discrepancy between the 
previous theories and the experiments, we theoretically 
investigate the quantum noise in LEDs at a low-injection 
level p2| . To this end, we extend the Langevin-equation 
method of Chow, Koch and Sargent |yj to treat the case 
where many photon modes are excited in an LED. From 
the microscopic quantum Langevin equations (QLEs) 
and the associated noise correlations, we derive the semi- 
conductor QLEs and the noise correlations at the LIL. 
An important assumption is that the photon-absorption 
and emission rates are much smaller than the photon- 
escape rate from a "cavity" . [] In the experiment |Q , it 
has been reported that rj^/rj > 2 in low-injection regions, 
whereas rj ~ rjd in high-injection regions. The differ- 
ence between rj and rjd is, in our theory, attributed to 
non-radiative processes and carrier-number dependence 
of lifetimes. They always exist in real LEDs, and be- 
come particularly important in low-injection regions. In 
contrast, we can easily show that the previous simplified 
theories always give r\ = r\&. Hence our theory is a 
minimal one to simulate real LEDs. A new formula for 
the photon Fano factor Q is derived. It gives the photon- 
number statistics as a function of the pump statistics, 
measuring frequency, rj, rjd and some factors arising from 
multimodeness of LEDs. Simplified formulae are derived 
for homogeneous cases, in which photons in each mode 
of the cavity are emitted and detected homogeneously, 
and for inhomogeneous cases, in which they are emitted 
and/or detected inhomogeneously. Using these formu- 
lae, we discuss the condition for generation of SPL in an 
LED. Our results agree with the experimental results. 

The paper is organized as follows. In Section 2, we 
derive the semiconductor QLEs and the noise correlations 
at the LIL. In Section 3, the semiconductor QLEs are 



The validity of this assumption with the meaning of the 
"cavity" will be explained in appendix A. 
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used to derive a new formula for the photon Fano factor. 
In Section 4, we examine the formula in various cases, 
and compare our theory with recent experiments [p|~[l2[ . 
In Section 5, we summarize the paper. 



II. DERIVATION OF QUANTUM LANGEVIN 
EQUATIONS AND THE NOISE CORRELATIONS 
AT A LOW-INJECTION LEVEL 

A. Microscopic Langevin equations of an LED 

Chow et al. j2l]] discussed the case where a single mode 
is excited among many-photon modes of a cavity. Since 
we are mainly interested in the photon statistics of LEDs 
at the LIL, we extend their method to treat many modes 
of photons. The total Hamiltonian 7i to t (which describes 
multimode photons in the cavity and carriers in the active 
layer of an LED) is written as 
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where 7i m uiti-ph is the Hamiltonian of the multimode 
photons, ai is the annihilation operator for the photons 
in mode I, and v\ is the field oscillation frequency in mode 
I. The Hamiltonian of the electrons and holes in the ac- 
tive layer is 7i C arricr- The annihilation operators of the 
electron and hole of wavevector k are Ck and <i_k, re- 
spectively, m c and are the electron and hole effective 
masses, respectively, and e° is the bare band-gap energy. 
The interaction among the carriers and the photons is 
represented by Hdipoic in the dipole approximation, with 
5/°k being the bare coupling constants, and h.c. means 
Hermite conjugate. Note that we here consider only a di- 
rect radiative transition, however, in real devices, there 
are other processes, and we will include them when nec- 
essary. The many-body interaction between the carriers 
is represented by 7i m any-body , the Hamiltonian of baths 
(or environments) by 7ibaths , and the interaction between 
the baths and the system (the carriers and photons) by 

^baths — sys ■ 

To obtain QLEs, we start with making a mean-field 
approximation for H m any-body, and consequently e° and 
g£ i ar e renormalized (See, e.g., Chap. 4 of Ref. pi[). 
The renormalized parameters are denoted by e g and gi ^. 

We then eliminate Wbaths + Hbaths-sys by using the 
Markov approximation |p|,|l6[-pT[| . As a result, the fluc- 
tuation and dissipation terms appear in the equations of 
motion. The generalized Einstein relation Jl6| |l8| , p0| , pl| 



gives the relation between the fluctuation and dissipation 
terms as follows, 



2D 



(5) 



where 2D^ V is a diffusion coefficient, a system vari- 
able, D M a dissipation term in Langevin equations, i.e., 
= D^+F^ (F M (t)F v (tf)) = 2D tlv 8{t-t l ). The brack- 
ets mean the ensemble average for fluctuations. 

We finally obtain the following microscopic QLEs, 
which describe LEDs in a microscopic scale, for the dipole 
operator Ok = d_kC k e Wi *, for the electric field operator 
Ai(t) — ai(t)e Wlt , and for the electron occupation prob- 
ability in k space n ck = c^Ck, 



dt 



Ok = -(7 + iwk - ivi)(Tk 

+»y^ 91' ,*Ai> (n e k + "hk - 1) + F ak , 



(6) 



d -° 



dt 
d 



k 

"ck 

5^* k 4a k +h.c.) + F ek , (8) 



,.-™ek = -Pek(l - "ck) - 

at t, 



where 7 is the dipole decay (dephasing) rate, fko^ 
1 .. + h 2 k 2 /2m e + h 2 k 2 /2mi l is the transition energy, and 
F^k is the fluctuation operator for the dipole. The hole 
occupation probability in k space is n^k = ^Lk^-k- The 
photon escape rate from the cavity is = vi/Qi where 
Qi is the Q factor of the cavity, is the passive-cavity 
frequency, and Fi is the fluctuation operator for the elec- 
tric field. The pump rate due to a current injection or 
optical pumping is F c k(l — n ck) where the factor (1 — n ek ) 
represents the pump blocking H]. Note that a lifetime 
of non-radiative decay r nr has been introduced in Eq. 
(H). As discussed later, the existence of non-radiative 
processes is, in our model, a necessary condition for the 
difference between the quantum efficiency and the dif- 
ferential quantum efficiency to occur. The other condi- 
tion is that the lifetime of radiative processes or that of 
non-radiative ones varies with n c . [See Eqs. (p6|)-(|68|).] 
Non-radiative processes might be modelled by capture of 
carriers at a trapping level in an LED, however, we only 
need that r nr is an implicit function of n c . These terms 
of pump and non-radiative decay have been phenomeno- 
logically introduced. The fluctuation operator for the 
electron number is F ck . 



B. Adiabatic approximation 

To derive more useful forms for later discussion, we 
use the adiabatic approximation |l^,^3|], and approxi- 
mate the solution of Eq. (H) by 
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i Hv 9i'M A v (™ek + «-hk - 1) + F ak 
0"k ^ — ; , (9) 

7 + tujk - ivi 

Substituting Eq. (|) into @), we find 

Ai = -[ K ^/2 + i(Sli-v,)]Ai 

+ J2Gii' A i'+ F i + F *,h (10) 
v 

where Gw is a "gain matrix" ; 

G w = ^5i*kS , i',kI'/,k(nek + ^hk - 1), (11) 

k 

where 2\k is a complex Lorentzian; T>i k 



7+i(wk— vi) ' 

A new fluctuation operator F a j (t) has been defined by 



(12) 



which is associated with the coupling between the carriers 
and photons. 

From Eq. (jlo|) , we also find the QLE for the photon- 
number operator rij = a\Ai, 



d_ 

dt 



n, = -K°ni + Y^[G IV A\A V + h.c] 



[(F* +F/)^+h.c], 



(13) 



Substituting Eq. (||) into @ , we also find the QLE for 
the total electron number n c = ^ k n e k> 



-^"c = 51 f e k(l - «ck) - 



2[G U /i4 I t i4 l /+h.c] 

i,v 

^^-^[il^+H, (14) 



C. Microscopic noise correlations of an LED 

To discuss the statistical properties of light emitted 
from LEDs, we must determine the noise correlations. 
Hereafter we assume that the correlations between dif- 
ferent modes of the photons and those between different 
wavenumbers of carriers can be neglected |19|, i.e., 



{A\(t)A v {t))~{m)5 w , (15) 

(0k(*)°v(*)) - ( n ck«hk)^k,k', (16) 

(<7k(i)4'(*)> - ((1 - "ck)(l - n h k))«k,k', (17) 

(nek(t)n e u.'(t)) - (^ek)<5k,k'- (18) 

We can then calculate the noise correlations, which are 
consistent with the QLEs @-(§), using (||); 



(F^(t)F ak ,(t')) ~ 2 7 (n ck nhk><5k,k^(i - t'), (19) 
(^k(t)4(0) - 2 7 <(1 ~ « ck )(l - n hk )) 

<Jk,k^(* - 1'), (20) 

(J^(t)F,/(0) =s itifHyfiSivSQ-H), (21) 

(fl(*)^(0) - K ? + ~ 0. ( 22 ) 

(F ek (t)F ek ,(t')) = [<P ek (l - u ek )) + (n ek /r nr )] 

<5k,k'«(t - 0. (23) 

where n{yi) is the number of thermal photons in mode I 
in the cavity. We have used, in Eqs. (|l^) and (^0|), the 
quasiequilibrium conditions [pl| , 



27<«ck«hk) > ^("ck^hk), 



(24) 



27((1 - «ek)(l - ™hk)) > -<(1 - n ek )(l - n hk )). (25) 

at 

We can also calculate the correlations of F CTi z; 

(Fl l (t)F a j(t , )) = (B^ l )S(t-t') 1 (26) 

(27) 



(^,j(*)^l(* / )) = (^ata,j)*(*-0» 



where R S p,i (Rabs,i) denotes the spontaneous emission 
rate into photon mode / (the absorption rate of photon 
mode I); 



Rsp.l = — = — ^ |ff;,k| 2 ^/,k«ek«hk, 

t t ,i 7 



(28) 



^abs,z = - l5i.k| 2 £/. k (l - n e k)(l - «hk)- (29) 



We have used Eqs. (19) and (|20|), and defined t t j as 
a radiative lifetime of the carriers into mode /. The 
(dimensionless) Lorentzian line-shape function £/k is 



C 



At this stage, we have the following QLEs, 
jni = -mm + [{Fl tl + F})Ai + h.c], 



(30) 



37 n c = ^ Pek(l - «ok) - — 
at T "— 



(31) 



k i 

where we have used 

Gil + G*; = i? SPi ; — -Rabs,Z, (32) 

and defined the renormalized photon escape rate 

k-i = — R sp ,i + -R a bs,i- (33) 
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D. Semiconductor Langevin equations and the noise 
correlations for an LED at a low-injection level 



In what follows, we assume ^ 

Rsp,l, Rabs,l <C Kl ■ 



(34) 



From Eqs. (]3G)) and (gj), the QLE for m at the LIL 
becomes 



a q 

—m = -n l m + R sp j + F n j , 



(35) 



where we have defined a new fluctuation operator; 
F n>l = [(fj , + F/)A ; + h.c] - («?n(^) + i? sp ,;)- (36) 

The right-hand side is averaged to be zero because the 
noises, F ni i and F& t i, have the Markovian property: 
(F?A t + h.c.) = Kfn{m) and (F^Ai + h.c.) = R sp>l . 

From Eqs. (|34]) and (|35|), the steady-state value of the 
photon number becomes 



(nj) 8 . s . ~ Rspj/tf < 1, 



(37) 



i.e., the photon number in each mode is quite small. 

On the other hand, from Eqs. @, @, and @, the 
QLE for n c at the LIL becomes 



d „ n c n c „ 
-n c = P- — - — +F c 

CUV Tp Tjjp 



(38) 



where P = ^ k -P e k(l — ™ek) is the total pump rate, and 

TV 71 _ fa ^ T i ' 



(39) 



is the radiative decay rate (r r is the radiative lifetime of 
carriers). The fluctuation operator for the total electron 
number F c is denoted by, 

f c = [-(4 F *,i + h - c -) + ^p,i] ( 4 °) 



Note that the lifetimes (t t ,i and r r ) are carrier-number 
n c dependent because R S p,i is an implicit function of n c . 

We thus obtain the final forms of the semiconductor 
QLEs at the LIL as follows, 



d n c n c 

— n c = p h r P + r r + r nr , 

dt r r r nr 



— ni = -rfnt + — + F K i + F r ,i, 

Ut Trl 



(41) 
(42) 



2 The following inequality, at first sight, might look rather 
strange. We therefore explain its validity at the LIL in the 
appendix A. 



where the noise operators, F n i,F c , have been divided 
into new ones, r r , F Tj i, r nr , Tp, F K< i. They are given by 



r r E^[-(4F ff ,, + h.c.) 



R 



•sp,i 



F Tj = A}F a ,i + h.c. - R 



sp.Z i 



K,l 



a\f 1 + f}a u 



(43) 

(44) 
(45) 

(46) 



where we have neglected the number of thermal photons, 
n(y{), since this is the case in usual experiments. These 
noise operators (r r , F T j, T nT , Tp,F K ,j) are associated with 
the radiative decay of carriers, the conversion from car- 
riers to photons, the non-radiative decay of carriers, the 
(intrinsic) pump fluctuation, and the photon escape from 
the cavity, respectively. The mean values of the noises 
are defined to be zero, of course. 

From above the noise correlation for T r can be calcu- 
lated as, 

(Tr(f)rr(f)) 

i 

+ (A c MAl l {t'))(Fl l (t)F a At')) 
= ^[(i? sPl ; + Rab a ,i){ni) + (R sp ,i)]S(t - t') 



J2(R S pMt~t') 



5(t-t'), 



(47) 



where, as done in Eq. (|l5|), we have neglected, in the 
second and third lines, the correlations between different 
modes. We also have assumed the following, 



r r ~^[-(At^ + h.c.) 



where A\ l is an operator uncorrelated with F a j [ pi) . 
Thus the property of fluctuations, (r r ) = 0, still holds. 
Equations (|26|) and ( p7j ) have also been used in the fourth 
line of Eq. Ml), and Eqs. and © in the fifth 

line of Eq. (47). 



Similarly we obtain 

(F Itl (t)F Ttl ,(t')) = 

(r nr (t)r nr (t')) = 



-)s(t-t'), 



— )5 u ,6(t~t'), 



6(t-t>), 



(T P (t)T P (t')) = (P)S(t-t'), 



(48) 

(49) 

(50) 
(51) 



where we have divided the correlation (r nr (i)r nr (i')) + 
(rp(i)rp(i')) into two parts; one is the correlation of the 
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non-radiative processes, and the other that of the pump. 
These semiconductor QLEs, ( |4l| ) and (^), and the noise 
correlations werejust assumed by us |2^,^5|] and 

Yamanishi and Lee (YL) H , although YL did not include 
non-radiative processes. 

Furthermore, we focus on the total photon flux, N, 
which can be detected at the PD surface. To this end, 
we must find two relations: one is the relation between 
the photon number in the cavity and the photon flux from 
the cavity. The other is the relation between the photon 
flux from the cavity and the photon flux detected at the 
PD surface. The input-output formalism ||,|2(J gives the 
first relation as, 



Vi 



F K 



(52) 



where Vi is the photon flux of mode / from the cavity. 
We next model the second relation using the argument 
of a beam splitter [^0[ . The relation between the mean 
flux of N and that of Vi is given by 



(53) 



where £/ is a "transmission coefficient" of mode / (see 
Fig. f), 

FIG. 1. Definition of the . Photons of mode / are emitted 
into the region between the solid lines. The photodetector 
can detect only some region (denoted by the one between the 
dashed lines) of the photons because its surface area is limited 
and the quantum efficiency rj of the detector is less than one. 
£z is the probability that a photon of mode I is detected by 
the detector. It increases as the surface area and/or rj do. 

Po is the transfer efficiency [||, and V = ^ z V; is the 
total photon flux from the cavity. The relation between 
the flux correlation of TV and that of VJ is given by 

(|AY(0)| 2 )=^(i-6)<Vi(tt)) 

+ £)6&,<AvUn)Afy(n)>, (54) 

where JV(f2) and VJ(f2) denote Fourier components of N 
and Vi, respectively. 



III. CALCULATION OF THE PHOTON FANO 
FACTOR 

In this section, we calculate the Fano factor of pho- 
tons H which denotes the normalized fluctuation of the 
photon number detected at the PD surface. 

Following the standard small-signal analysis used in 
Refs. [|]|],|7],[l4|] , we expand the radiative and non- 
radiative lifetimes, T r j[n c ], r r [n c ], and r nr [n c ], to linear 
order in An c = n c — n c o] 



Tr,i [n c \ 
r r [n c ] 

T n r[n c ] 

I 

TrO 



(r riI )o(l-«r,I— ), 
An c 

T nr o(f + K m ), 



n c0 



7"r0 



(55) 
(56) 
(57) 
(58) 
(59) 



where quantities with the subscript denote the average 
values at a running point, P — P . Sensitivity of lifetimes 
to the carrier-number fluctuations An c is represented by 
K r i, K r and K m . Though K t is usually positive, it is not 
always the case as we will discuss in the next section. 

Linearizing Eqs. (41), ( [42| ) and ( fffi ) in terms of An c , 
Am =m- (n/) 0) AP = P - P , and AV^ = Vj — (VJ) , 
and substituting Eqs. rt5q) and (p7\) into them, we obtain 



— An c = AP 
dt 



r r 



d . n . An r „ „ 

—Am = -k Ah, + — + F K i + F r ,j, 
at tL , 



(60) 

(61) 
(62) 



AVi = K^Ani - F K , U 
where we have used the following equilibrium conditions 

/ 1 1 X 

Po = ?^co( 1 ), 

T tO T aiO 

En, s n c0 n c0 

«"w)o = V t — r = — 



J2 v m = v ° 



Po 



1 + (r r o/r nr o) 
and introduced the effective lifetimes defined as 

/ r nrO 



(63) 

(64) 
(65) 



Dropping the noise terms in Eqs. (|60|)-(|62j), and using 
Eq. (|5^), we calculate the quantum efficiency rj and the 
differential quantum efficiency ?7d as 



TrO 




Ki)o 


1 + K t ' 


1 + K It i 




I 


I I 
~ + ~ 






r nr 



1-K n 



Po/t, 



r0 



00 



No = 

Po (1/Tro + l/r nr0 ) l + e ' 



where 



AN(n) 



AP(Q) 



Po 



(l/r r ' + 1K T ) 



TrO 
TnrO 



1-K m 

1+K r 



e . 



(66) 
(67) 

(68) 
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These efficiencies (rj and n^) are illustrated in Fig. 2. 

FIG. 2. Typical I-L (Injection-Light) characteristics of 
LEDs are shown by the solid lines. The quantum efficiency 
n is the slope of the straight line (represented by the dashed 
line) connecting the origin and the running point, whereas the 
differential quantum efficiency r/d is the slope of the tangential 
line (chain line) at the running point. 

Note that the difference between rj and fyj, which was 
measured to be large (e.g., r/ d /r/ > 2 at the LIL can 
not be explained by the previous theories which assumed 
T nr o = oo || ^| or K r — K nl = 00. In this sense, our 
model is a minimal one that simulates real LEDs. 

Hereafter we assume that il <C (~ 10 12 Hz), which 
is well satisfied in usual experimental conditions. Then 
the response of the output photon flux to the modulation 
of the pump is calculated as, 



AN(n) 



AP(O) 



(69) 



Hence 1/r" is understood as a cutoff frequency. It is 
seen that the non-radiative processes push the cutoff to 
the higher frequency [See Eq. (pq)]. 

On the other hand, the photon flux correlation between 
mode I and m is calculated using Eqs. (|60|)-(|6^) as 

(AVr(0)AF ro (0)) 



1 + (Sir") 2 

•' .„ i + in T " 



(70) 



where a total pump noise AP to t has been defined as 

AAot(fi) = AP(fi)+f P (fi), 

which consists of the modulation AP(il) and the intrinsic 
pump noise Tp (O) . Note that the latter noise can be sup- 
pressed by inelastic scattering in conductors ]4|J^- [l0| , |26|| . 

To see the physical meaning of Eq. (70), let us intro- 
duce the Fano factor of the pump electrons (or excitons) 
W e and that of the p/iotons detected at the PD surface 
Wph, which are defined by 



W e (Q) 



(|AP to t(»)| ; 

P T 
(|AiV(ft)| 2 ) 
N T 



(71) 
(72) 



where T is the Fourier-integral time. We transform Eqs. 
mT, ° the Fourier components: 

(P;(n)P r (n)) = <f r *(n)f r (fi)) = — <f * r (Q)f nr (^)> 

£0 

= -(f?(0)F r (0)) = — T = V T= (73) 
T r0 1 + e 



Substituting Eqs. (|7(]) and (|7|) into (||), and dividing 
by NqT, we finally obtain 



w P h(fi) = i - J 



4 1 + W e (Q) 



where we have used Eqs. (p3[), (|63[), ( p6| 
the following, 

C — \ " r r Tr0 _ 
f^r r ',(r r>ro )o/3o/3o 



+ (ftr") 2 ' 77 l + (f2r") 2 ^' ^ 



381) , and defined 



(75) 
(76) 



which represent effects due to multimodeness of LEDs, 
and the consequences are discussed below. Note also that 

Ci < i,C 2 < i- 

Equation y74j) is our main result which gives the Fano 
factor of the photons detected at the PD surface as a 
function of the Fano factor of the pump, measuring fre- 
quency, and several parameters (r\, rjd, £i, C2, and r"). 
Note that 1 jr" becomes a cutoff frequency of the Fano 
factor as well as that of the modulation [see Eq. (|69|)]. 



IV. DISCUSSION 

A. Low-frequency limit 

Let us examine Eq. ( ]74| ) in two cases, VL — and £1 > 0, 
separately. We first discuss the low- frequency limit. This 
case is applicable to the most experiments, because they 
are usually performed at the frequency which is lower 
than any other relevant frequencies. Setting Q — * 0, we 
obtain 

VFph(0) = 1 - 277^1 + —[1 + W c (0)}( 2 . (77) 
V 

There are several cases where this expression becomes 
simpler: 



1. The case where photons are emitted and detected 
homogeneously (the homogeneous case) 

(h-1): When photons in each mode are emitted and de- 
tected homogeneously, i.e., when K r ^i = K x — constant 
and = /3q = constant, we have £1 = £2 = 1, hence 



Wp h (0) = 1 - 2? 7d 



V 



W e (0)}. 



(78) 



This is the formula which we have derived elsewhere p5[ . 
We have implicitly assumed this property (h-1) there, 
and this is a natural choice unless one considers a situ- 
ation such that inhomogeneity due to, e.g., cavity-QED 
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effects becomes important. 

(h-2): In addition to (h-1), when the non-radiative pro- 
cesses do not exist and/or the carrier-number dependence 
of lifetimes cancels to be zero, i.e., r nr o — > oo and/or 
K T + K nT = 0, we have 77 — rjd and the I-L (Injection- 
Light) characteristics become straight. Hence, from Eq. 
(|7q), we obtain 



W ph {0) = 1-2t? 

= 1-TI + 



hr?[l + We(0)] 

»?We(0). 



(79) 



This is just the previous formula which is frequently used 
in the literature §-|@. 



2. The case where photons are emitted and/ or detected 
inhomogeneously (the inhomogeneous case) 

Note that in the homogeneous case the factors £1 and 
£2 do not appear in the expressions for W p h [Eqs. ( |7g| ) 
and (|79|)1. On the other hand, they appear when emission 
and/or detection efficiencies are different among different 
modes. We call this general case the "inhomogeneous 
case." To investigate this case, let us consider a simple 
case where K T j = and K nl 7^ 0. In this case, (1 — (2 = 
£, where 



E 



TlQ & 
(Ti,i)o Po 



Hence, from Eq. 



we obtain 



W ph (Q) = 1 - 2 Vd C + -£Si[l + Wc(0)] 



1 



(80) 



(81) 



Compared with Eq. (|78|), we see that 77 and ??d are ef- 
fectively multiplied by £ in Eq. (81). However, this does 
not mean that C could be absorbed in r\ and 7yd, because 
they are already defined by Eqs. (|6^) and (|67|), respec- 
tively. Any redefinition would lead to disagreement with 
the observed I-L characteristics. 



3. Condition for generation of sub-Poissonian light 

As another illustration of our result, we next dis- 
cuss the condition for generation of sub-Poissonian light 
(SPL) (W p h < 1) with a Poissonian pump (W c = 1). For 
simplicity, we hereafter consider the case (h-1) [Eq. ((78])], 
because many LEDs seem to be categorized in this case. 
From Eq. (J7§), we obtain 



< 77 d < 77, 



(82) 



as the condition. This can be intuitively understood that 
the flatter the I-L characteristics, the duller the sensitiv- 
ity of the LED to the pump fluctuation. Hence, if I-L 
characteristics are like Fig. 2 (b), then even a Poissonian 



pump (W e = 1) can produce sub-Poissonian light. Note 
that this mechanism is completely different from those of 
Refs. @,[l|,|7|, because the authors of |@|1||7) eventu- 
ally make the current noise injected to the active layers 
below Poissonian. On the other hand, the condition ( |82] ) 
is equivalent to 



K T + K BI < 0, 



(83) 



where we have used Eqs. (p6[)-(p8|). For simplicity, we 
here take K nr = 0, and find below what is needed for 
K t < 0. 

The spontaneous emission (SE) rate n c /r r can be ex- 
pressed approximately as n c /r r cx (n c ) p , where p is a 
constant. We therefore have K T ~ p — 1 from Eq. (|5^). 
It is well known that p ~ 1 (p ~ 2) in high- (low-) 
injection regions for SE processes of free carriers. For 
exciton recombination, we have p ~ 1 at the LIL. Thus 
we usually have non- negative values of K T . Making use 
of cavity-QED effects, however, we can obtain negative 
values of K r . This is illustrated in Fig. 3 for a p-doped 
quantum well structure in a micro cavity, where we have 
assumed the following: the conduction band is parabolic, 
the effective electron mass is 0.1 times the free electron 
mass, the doping level is high, and the cavity-QED effects 
prohibit SE except at the band edge. 

FIG. 3. The spontaneous emission (SE) rate as a function 
of carrier population in a doped-quantum-well structure in a 
micro cavity. Three curves correspond to different tempera- 
tures: T = 80K (bottom), T = 15K (middle), and T = 3K 
(top). 

It is seen that K T becomes negative for the sheet carrier 
density > 10 14 m -2 when temperatures are low enough 
(maybe below 3K). In this case, we obtain r]d < rj [see 
Eqs. (|66[)-(|68|)], and even a super-Poissonian pump can 
produce SPL. Exciton recombination with cavity-QED 
effects would work better, which will be discussed else- 
where. 

It has been therefore shown that our formula ( |74| ) is 
useful to find the condition for generation of SPL, or, of 
course, other quantum states of light. 



B. Finite- frequency cases 

Let us turn to the finite frequency cases. For simplicity, 
we here use Eq. ( |T4| ) with property (h-1), i.e., Ci = C2 = 
1, and assume the case where K T = K nr in Figs. 4. 
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FIG. 4. The photon Fano factors are plotted as functions 
of frequency Q (in units of I/tvo). (a) and (b) represent the 
cases where K T = = 0.5. (c) and (d) the cases where 
K T — K nT — —0.5. The pump noise level is assumed to be 
zero (We = 0) and the Poissonian level (W e = 1) in (a), (c) 
and (b), (d), respectively. In each figure, the solid and dashed 
lines represent the cases where the non-radiative recombina- 
tion is significant (r r o/r nr o = 1) and absent (r r o/r nr o = 0), 
respectively. The cutoff frequencies are indicated by vertical 
arrows. 

The solid lines (dashed lines) in Figs. 4 represent the 
cases when non-radiative processes exist (do not exist). 



1. When K T = K nr = 0.5 

The noiseless (W c = 0) and Poissonian (W c — 1) pump 
are denoted by Figs. 4 (a) and (b), respectively. When 
T nr o = oo (dashed lines), we recover the results similar 
to those of Yamanishi and Lee § . We also see that the 
cutoff frequency (indicated by vertical arrow) for r nr o 7^ 
00 becomes higher than that for r nr o = 00 [Fig. 4 (a)]. 



2. When K T = K n 



-0.5 



The noiseless (W e = 0) and Poissonian (W e = 1) pump 
are denoted by Figs. 4 (c) and (d), respectively. We find 
that the existence of non-radiative processes gives smaller 
Wph in some frequency region [Fig. 4 (c)] or in the entire 
frequency range [Fig. 4 (d)]. In particular, Fig. 4 (d) 
shows that the Poissonian pumping can produce SPL in 
a wide frequency range. Figure 4 (c) also shows that the 
cutoff frequency for r nr o ^ 00 becomes higher than that 
for r nr o = 00 as well as Fig. 4 (a) does. 

It might be difficult, however, to observe these fea- 
tures experimentally, because they would not appear up 
to l/T r o ~ 1GHz. (Recently the photon Fano factor of 
LEDs has been measured up to 40 MHz |12|.) 



C. Comparison with the experimental results 

Hirano and Kuga M reported that the previous for- 
mula (|7^) disagrees with their experimental data in low- 
frequency regions: They measured the following ratio 



W ph (W c = 0) 

Wph (We = 1) 



(84) 



and obtained that r < 1 — 77, whereas Eq. (|79| ) gives 
r = 1 — 77. Since the measuring frequency r2 mca s is low 
enough, i.e., f2 meas ~ 10MHz <C 1/r" ~ 1GHz, and their 
LEDs seem to be categorized in the homogeneous case 
(h-1), we can use Eq. (|78j ) to analyze their results. Our 
formula (78) gives 



1 - 2 m + rip 7] 
1 -2^ + 2772/77 ' 



(85) 



When ?7d > 77 (which is usually the case in low-injection 
regions [p|-^2[), our theory gives r < 1 — r\, in agreement 
with the experimental results. Furthermore, Hirano et 
at recently confirmed that Eq. p^) agrees with their 
experimental data (see the table below). It is seen that 
the values of the previous theory (r = 1 — 77) are larger 
than the experimental values, whereas the values of Eq. 
(|S5| ) are closer to them. It also has been shown that, 
when W c = 1, W p h itself is more than one [jllj. This fact 
also supports our formula ( f78| ) because it gives, when 
W e = 1 and 77 d > 77, W p h = 1 - 2?7d + 277^/77 > 1, and 
their LEDs also have the property 77d > 77 Pjllj . 

On the other hand, a simple argument flll|,p8| leads to 



Wph(0) = l-r7+^We(0). 

V 



(86) 



Since |?7 — 77a | is small in Table I, the difference between 
this formula and Eq. ( |78| ) is within the experimental error 
and is not detectable. Further experimental studies are 
needed to observe the difference. 



V. SUMMARY 

From the microscopic QLEs ([|)-(|]), the associated 
noise correlations (19)-(p3|), and the assumption (34), 
we have derived the effective semiconductor QLEs (41), 
( [42| ) and the associated noise correlations (f47|)-(|5"l|). The 
assumption (|34|) is valid at a low-injection level and in 
real devices as explained in appendix A. Applying the 
semiconductor QLEs to semiconductor LEDs, we obtain 
a new formula J74) for the Fano factor of photons. It 
gives the photon-number statistics as a function of the 
pump statistics and sever al p arameters of LEDs, which 
are defined by (|65|)-(|6^), ( |7E| ) and (|76l). Key ingredients 
are non-radiative processes, carrier-number dependence 
of lifetimes, and multimodeness of LEDs. The formula 
is applicable to the actual cases where the quantum ef- 
ficiency 77 differs from the differential quantum efficiency 
rid, whereas the previous theories ^-JtJ turn out to give 
V = Vd- It is also applicable to cases where photons in 
each mode of the cavity are emitted and/or detected in- 
homogeneously. When 7/d < 77 at the running point, in 
particular, our formula predicts that even a Poissonian 
pump can produce sub-Poissonian light (see Sec. IV A 3). 
This mechanism for generation of sub-Poissonian light is 
completely different from those of the previous theories, 
which assumed sub-Poissonian statistics for the current 
injected into the active layers of LEDs. It was shown 
that our results agree with recent experiments by Hirano, 
Kuga and coworkers [p|-|l2||. We have also found that, 
in finite frequency regions, non-radiative processes some- 
times give better results (smaller W p h and/or a higher 
cutoff frequency). These will deserve further theoretical 
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and experimental researches of quantum aspects of light 
emitted from LEDs. 
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FIG. 5. (a) A schematic band diagram of real LEDs at the 
LIL, (b) the Stokes shift [CB (VB) represents the conduction 
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"cavity". 
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APPENDIX A: VALIDITY OF THE INEQUALITY 

We here show that the inequality (|34|), i? sPj z, i? a bs,z <C 
Ki, holds for LEDs at a low-injection level (LIL). 

At the LIL, pumped carriers first relax to lower-energy 
states which are formed by impurities, defects, spatial 
randomness and so on, and then recombine to radiate 
photons [Fig. 5 (a)]. Within the energy region of our 
interest, where photons are emitted, we thus have R sp j ~ 
-Rabs,z- I n addition, in real devices, there exists a Stokes 
shift [Fig. 5 (b)] which makes, at low temperatures, i? a bs,z 
even smaller. Hence we have the absorption and emission 
profiles such as Fig. 5 (c). 

We now proceed to compare i? S p(abs),z with In Eqs. 
(H) and (p9|) ; |gz, k | 2 is proportional to 1/V cavi ty (Kavity 
is the volume of a cavity), and ^ k is to the volume 
of the active layer Votive, hence we have i? S p(abs),z oc 
^activc/^cavity I n the analysis of laser diodes (LDs), 
it is customary to take Votive — Cavity j because the 
devices are made so that the lasing modes are con- 
fined in the active layer. In this case, we usually have 
Rzbs.i/c ~ 10 4 cm _1 > rf/c ~ 10 2 cm _1 (c is the ve- 
locity of light), thus we can not obtain the inequality 
([m|). In the case of LEDs, on the other hand, they are 
usually designed in such a way that the reflection coef- 
ficients on the boundary surfaces are small, hence most 
modes of photons of our interest are not confined in the 
layer. In this case, it is natural to take the "cavity" 
volume as big as a cube on which the detector's sur- 
face is located [Fig. 5 (d)] . Then k° can be estimated as 
1/k° ~ (V r cav ity) 1/3 /c+Qtdcvico, where ^device is a time for 
photons to traverse the device. S ince Q and ^device are 
small for LEDs, we have k° ~ c/ (Verity) 1 ^ 3 ■ Therefore, 
noting that i? abs ,z oc 1/T4avity and k° oc l/(14 avit y) 1/3 , 
we can have the relation i? a bs,z *C k®, if we take V^ av ity 
big enough. Thus we obtain the inequality ([mJ) for the 
case of LEDs at the LIL. 
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TABLE I. The experimental values of r),r]^,r and the the- 
oretical values of r. After Ref. [11]. 



?7(exp.) 


r?d(exp.) 


r(exp.) 


r(Eq. (35)) 


r = 1 — r\ 


0.067 


0.090 


0.90 


0.89 


0.93 


0.104 


0.125 


0.84 


0.86 


0.90 


0.150 


0.175 


0.81 


0.81 


0.85 
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